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Abstract 

We show that threshold resummation of the end-point logarithms in inclusive 6-hadron decays leads to a Sudakov 
factor, which introduces l/rrifc power corrections through infrared renormalons. These power corrections are strong 
enough to account for the observed semileptonic branching ratio, charm yield, and lifetime ratios by breaking the 
local quark-hadron duality. 

Several puzzles in inclusive heavy hadron decays have existed for some time. The observed semileptonic branching 
ratio BgL and charm yield Uc of the inclusive B meson decays were Bsl — (10.19 ± 0.37)% and ric = 1.12 ± 0.05 
from CLEO §, and Bsl = (11-12 ± 0.20)% and Uc = 1.20 ± 0.07 from LEP g. For updated values, refer to 
1^. The naive parton model, which coincides with leading-power heavy quark effective theory (HQET) gives 
Bsl - 13% [|. The 1/m^ corrections are less than 5% of the leading ones M. Another puzzle is the low lifetime 
ratio r(Ab)/r(Brf) 0.794 ± 0.053 g. The HQET prediction up to 0{l/ml) is about 0.99 §. When including 
the 0(1 /ml) corrections, the ratio depends on six unknown parameters, and reduces only to around 0.95 for various 
model estimates g. 

There has been the suggestion that l/rrifc corrections might appear, i.e., the local duality might break down in 
nonleptonic decays It has been observed that if rrih was replaced by the 6-hadon mass in the phase-space factor 
TO^ associated with nonleptonic decays, the predictions agree well with the above data This modification reduces 
Bsl by increasing the nonleptonic branching ratios, and explains the absolute B meson decay rate [0. Note that 
HQET based on the b quark kinematics accounts for only 80% of the decay rate. However, the mass replacement is 
lack of a solid theoretical base. 

In this letter we shall propose possible mechanism responsible for the breakdown of the local duality, which arises 
from threshold resummation associated with final-state particles. It has been shown that resummation for the intial 
heavy hadron does not violate the duality [ll| ] . Define the variable A — niB — nih as the B meson and b quark mass 
difference. In the kinematic region with the invariant mass of the outgoing quark being O(Amh), infrared divergences 
in the decays are factorized into a B meson structure function /(fc) ]l2| , p^ , k being the residual momentum carried by 
the b quark. The first moment of f{k) vanishes, implying that nonperturbative corrections start from 0(l/m^). In the 
end-point region with the invariant mass vanishing like O(A^), additional coUinear divergences from loop momenta 
parallel to the energetic quark momentum occur in higher-order corrections, which demand the introduction of a jet 
function. Take the b —>■ sj transition as an example. The end-point logarithms asln(l — x)/{l — x)+ (as In^ iV in 
Mellin space) with x — 2E-y/mf), where E-y is the photon energy, should be orgainzed into a jet function J(l — x). 

Threshold resummation of the above double logarithms leads to |14] 
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where c is an arbitrary real constant larger than the real parts of all the poles in the Sudakov factor J{N), and the 
anomalous dimension is given, up to one loop, by "/k = cusCfIt^, Cf — 4/3 being a color factor. Note that the 
physical range of a; is < a; < mB/iTii,. We show that infrared renormalons in Eq. (Q) generate l/rri), corrections by 
inserting the identity as{n^) = tt J^^ da exp[—aj3i \n{iM^ / Aqqj^)] into J(Af), with the coefficient /3i = (33 — 2ri/)/12, 
Uf being the flavor number. Performing the integration over A, we have 
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Equation (g) can be expressed as a sum of contributions from all the infrared renormalons. It is known that 
the leading infrared renormalon from a —f corresponds to perturbative resummation. This can be understood by 
considering the a ^ limit of the integrand in Eq. (0) : 
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The integrand on the left-hand side has the same functional form as in J{N), if as was fixed, except the undetermined 
const.. Hence, the left-hand side is equivalent to J{N)\a^ at the leading- logarithm accuracy. Substitute Eq. (||) into 
the inverse Mellin transformation in Eq. (|l|), and choose the branch cut along the negative real axis in the complex 
iV plane. For < a; < 1, we enclose the contour in Eq. (|l|) toward the minus real axis. The discontinuity crossing the 
branch cut gives, using the variable change w = \nN [n5f. 
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because the integrand is an odd function of w. For x > 1, we enclose the contour toward the positive real axis, where 
there is no cut and no pole, indicating J{1 — x)\a^ = for x > 1. It is obvious from Eq. (|l|) that J(l — a;)!^^ vanishes 
at a; = 0. Therefore, we obtain the qualitative behavior of J(l — x)\a^ : it takes a value for < a; < 1, and diminishes 
at a; = and at a; = 1. 

The next-to-leading pole from a l/(2/3i) leads to 
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We neglect the second term in the brackets, which is not important compared to the first term, and employ the 



approximation z 



(1 + z-l) 



Af-1 



1 + N{z — 1), since the integral over z is dominated by the region z ^ 1. 



Equation (|5|) becomes exp[const.(AQCD/™6)-^] ~ 1 + const. (Aqcd/™&)-^, up to corrections of ©(AqQ^/m^), where 
the divergent integral J dz / {1 — z) has been absorbed into the ambigious const.. Picking up the first two renormalons, 
we write 
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Inserting the above expression into Eq. (|l|), the second term leads to const. (AQCD/™6)a:^>/'(l — x)\as- To have 
meaningful predicitons, a 1/777;, power correction must appear to cancel this ambigious term. Viewing its functional 
form, this l/rrib correction can be introduced by Fourier expanding the jet function: 
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where the variable change x — {niB /nib)t, < i < 1, has been made. It can be shown that the a l/(2/3i) 
renormalon in the first term on the right-hand side of Eq. (j^) is cancelled by the cr ^ renormalon in the second 
term under an appropriate prescription for extracting the renormalon contribution. The details will be published 
elsewhere. After this cancellation, we freeze the coupling constants in the above two terms, and obtain well-defined 
predictions for inclusive 6-hadron decays. 

We study the semileptonic decays b — > civ for the lepton I — e,^,T. Adopt the scaling variables x — 2Ei/niB and 
y = q^/ni^ for the b quark decays, where Ei is the lepton energy and q = pi + p^, the momentum of the lepton pair. 
The physical bounds in the ranges. 
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can be achieved, since the b quark carries residual momentum [l^ . The ratios are defined by ai^D,c = TOi,n,c/™B, 
mi, niD and nic being the lepton mass, the D meson mass, and the charm quark mass, respectively, niu appears as 
the minimal invariant mass of the decay product Xc- 

For these modes, loop corrections do not generate double logarithms, because the c quark is massive and coUinear 
divergences are absent. The factorization formula for the total decay width is then as simple as 
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^ mi I dEidq'H ( ^ ) = To ( ^ ) I dxdyH [--^x/^y] , (9) 

with To = 01,1^6^771^/(1677^) and the hard part H{x,y) = (1 + y — x — ml/ml){x — y — mf/ml). Equation (||) 
does not involve the structure function f{k), since the higher-power contributions it introduces start from 0{l/mf). 
Expanding the above expression up to A/ms, we obtain 



Fw J dxdy^^ ^1 + {1 + y - x - al){x - y - af) 



-— [{y + af){l -2x + 3y - Saj) - x{2y ~x~ 2al)] ] . (10) 

For the nonlcptonic decays b —^ ccs, the threshold resummation effect is associated with the energetic s quark, and 
both the c and c quarks are regarded as being on-shell. The definitions of the scaling variables are modified into 
x = 2Ec/mB, y = (pc +Pc)'^/m%, and j/o = 2(i?c + Ec)/mB, with Ec (i?e) and Pc (ps) being the c (c) quark energy 
and momentum, respectively. Their physical ranges are 

2ac< x<l 

(1 - x){x - y/x^ - 4a2) + 2al ^ ^ ^ {l-x){x + y^x^ - 4a2) + 2al 



2{l~x + al) - 2{l-x + al) 

xy- y/x"^ - Aal^/ y{y - Aal) < yo< I + y . (11) 

The factorization formula for the total decay width is written as 
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Tcs = To / dxdydyoc{mb)H x, — jy, yo ] J Vo H , (12) 

\mb J J \ nib ml nib J \ mb m^ J 

with H{x, y, yo) — (yo — x){x — y). The factor c(mh) is given by c{nib) — {Nc + l)c^(TOb)/2 + {Nc — l)c?_(mh)/2 with 
the color number Nc = 3 and c± — C2 ± ci, where the Wilson coefficients ci and C2 correspond to the four-fermion 
operators Oi = {sL"^^bL){cL^^CL) and O2 = {cLl^ibL){sLl^CL), respectively. We have confirmed that our predictions 
are insensitive to the choices of the arguments of the Wilson coefficients. Expanding H and J up to 0(A/mb), Eq. (|l^ ) 
becomes 

Tcs = ^0 y dxdydyoc{mb)^ + 6-^^ {yo - x){x - y)5{\ - yo + y) 

— — [{yo - x)yj{\ -yo+y) + {yo - x){x - y){yo - 2y)J'{\ - yo + y)] \ ■ (13) 
mb J 

Note that the jet function for the leading term, representing part of perturbative corrections, has been set to the 
(5-function. This is for consistency, because the semileptonic decays are evaluated to leading order in as- 

For the b ~^ cud mode, we associate the threshold resummation effect with the energetic u and d quarks, and regard 
the c quark as being on-shell. The scaling variables are defined hy w — p\/ni\, x — 2Ec/mB, y = {pc + PuY /m^^, 
and yo = 2{Ec + Eu)/mB, whose ranges are 

0< u;< {1-acf , 
2ckc < X < 1 ~ w + , 



X — J x^ ~ 4a'^ 2 ~ X + J x^ ~ 4a^ x + J x^ ~ 2 — x — J x"^ — 4a^ 
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The factorization formula for the total decay width is written, up to 0{K/mb), as 

Tud = ^0 J dwdxdydyoc{mb)l^ ^1 + {yo - x){'w + x-y- al)S{w)S{l - yo + y) 
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The prefactors To for the nonleptonic modes are the same as that for the b civ decays with \Vud\ = \ Vcs\ = 1- 

Setting vriB to rn;,, i.e., A = 0, Eqs. (^0|), and ([l5| ) reduce to the quark-level formalism. To estimate the 
strength of the A/m;, corrections, we propose the parametrization, 

Jp(l-t) = (a + l)(a + 2)r(l-t) , (16) 

where the free parameter a is determined by fitting Jp{N) to J{N)\a^ for the first few N . Equation ( p^ ) is normalized 
to unity, approaches the ^-function (5(1 — t), i.e., the lowest-order expression of the jet function, as a — > oo {as — > 0), 
and obeys the boundary conditions Jp(0) = Jp(l) = 0. For a fixed as, we have the moments J(2) = exp(— 7i<-), 
J(3) = exp(— 77i<-/4), • • •. The fit of Jp{N) to these moments implies a = 15 ^ 25. We have confirmed that the 
numerical results are insensitive to the variaton of a within the above range, and to the parametrization form of the 
jet function. Below we shall choose a = 20. 

The above formalism applies to the B^, Bg, and At, decays simply by redefining the variable A — mn — mi,. It is 
observed that the increase of the b quark mass just decreases the total decay rates, and has a small effect on the relative 
rates among different modes. The predictions are more sensitive to the variation of the charm mass: i?sL increases, 
while Tic decreases with rric. The total decay rates also decrease with nic. We choose mb = 4.8 GeV (correspondng 
to c{mi)) — 3.24) and = 1.5 GeV as representative parameters. The results from the quark-level calculation 
and from the inclusion of the l/nib corrections are displayed in Table I, for which the masses = ms^ = 5.279 
GeV, mB, = 5.369 GeV and mt^^ — 5.621 GeV have been adopted. It is found that the 1/mb corrections enhance 
the b — !■ cud decay rate more than the b — s- ccs decay rate, such that i?gL is reduced without increasing n^. The 
consistency with the observed lifetime ratios T{Bs)/T{Bd) and T{Kb)/T{Bd) is also improved. The ratio T{Bu)/T{Bd) 
remains to be unity, implying that the Bd and i?„ meson lifetime difference may be attributed to 1 /m^ or perturbative 
corrections |l^. Note that the absolute heavy hadron lifetimes for \Vcb\ = 0.04 are still larger than the data by 30%, 
thought the consistency is greatly improved. This discrepancy is expected to be removed by subleading contributions. 
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Table I. The predicted Bsl, n„ T{Bu)/T{Bd), T{Bs)/T{Bd), T{Ab)/T{Bd), and T{Bd)/T{BdUp for r(Bd)cxp = 
1.56 X lO^^g The data of the lifetime ratios are T{Bu)/T{Bd) = 1.066 ± 0.020 and T{Bs)/T{Bd) = 0.945 ± 0.039 |]. 
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